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One of the most dangerous parts of the flight is the landing phase, as most 
accidents occur at this stage. In order to reduce the effect of the low-level 
wind shear on the longitudinal motion of the aircraft in the glide path 
landing mode (task) a robust H— control isproposed. Dynamic models of the 
plane and wind shear are built. H, and Hæ synthesis methods are 
investigated for the task of aircraft flight control in a vertical plane during 
landing under conditions of undefined disturbances. Both control methods 
allow to reduce height deviation significantly. However, suboptimal control 
Hæ provides better quality of transition processes both in height and speed 
than optimal control H3. The results of simulation of the synthesized system 
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1. INTRODUCTION 

Ensuring flight safety is an urgent problem of modern aviation, especially during flights in difficult 
meteorological conditions. The most dangerous meteorological phenomenon for aviation flights is low-level 
wind shear with large gradients of wind components in height and range, which is caused by local 
disturbances of the atmosphere. Suddenly arising disturbances in the state of the atmosphere are extremely 
dangerous during an aircraft landing. They led, for example, to two known disasters: at New Orleans 
International Airport on July 9, 1982 the plane Boeing B-727 crashed during landing and at Dallas 
International Airport on August 2, 1985 plane Lockheed L-1011 crashed during landing. Due to the great 
relevance of the problem developers around the world were engaged in the construction of automated flight 
control systems capable of preventing such catastrophes. Various control algorithms based on different 
physical principles and mathematical concepts have been proposed, built for different models of the local 
state of the atmosphere, which somehow solved this problem [1]—[5]. 

The studies [1], [2], [5] consider the analysis and design of a robust controller. The controller is 
significant component of an entire automatic landing system developed as part of the aircraft landing task, 
which was proposed by AIRBUS and ONERA. Techniques of robust synthesis (e.g., structured Ho 
synthesis) acts as an effective basis for accomplishment of these tasks. In research [3] a robust automatic 
landing controller (SIRAC) based on stable inversion (SI) is proposed. The SI algorithm improves an 
indicator such as the output tracking accuracy, at the same time the application of the Hæ synthesis is aimed 
at increasing the robust stability to uncertainties that arise because of wind disturbances. In research [4] the 
vertical speed of the aircraft prior to landing is controlled based on a structured Hœ — control structure, 
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minimizing the effects of wind shear, ground effects and airspeed changes. A specific multi-model strategy is 
considered accounting changes in mass and center of gravity location. 

In research [6] a pseudo-sliding mode control synthesis procedure is considered and applied to 
develop control system for a nonlinear NASA Langley generic transport model. The synthesized control 
system allows minimizing aircraft loss-of-control by maintaining primary pilot input-system response 
characteristics throughout the flight, taking into account the possibility of actuator damage. In research [7] 
the problem of robust active fault-tolerant control (FTC) was considered for systems with undefined linear 
parameter variation (LPV) with simultaneous actuator and sensor failures. In research [8] a parameter 
independent embedded sliding mode controller with self-adaptation was developed that converges in the 
system in a finite time, focusing on the uncertain linear parameter variation (LPV) model of the aircraft 
variant, which has large-scale sweep angle variation and expansion. In research [9] a robust control scheme 
for commercial aircraft is presented. The control law is supplemented with a prior information about wind 
gusts. The main contribution of this article is the integration of the wind gust alleviation system using light 
detection and ranging (LIDAR) with the widely used flight control architecture the so-called C* control law. 
In research [10] the active vibration control of composite panels with uncertain parameters in a hypersonic 
flow is studied using the non-probabilistic reliability theory. Using the piezoelectric patches as active control 
actuators, dynamic equations of panel are determined by the finite element method and Hamilton’s principle. 
The results of the research prove the fact that the control method influenced by reliability, Hœ performance 
index, and approach velocity is effective for the vibration suppression of panel throughout an entire interval 
of uncertain parameters. 

The research [11] developed a multi-loop controller for a morphing aircraft to guarantee stability of 
the wing-shaped transition response. The suggested controller uses a set of inner-loop gains to ensure 
stability using classical techniques, whereas a gain of self-adaptive Hœ outer-loop controller is designed to 
provide a certain level of robust stability and performance of the time-varying dynamics. The paper [12] 
describes an analysis method, a generalization of the developed parameters of the robust controller for 
aircraft lateral control using auxiliary damping automatic devices (ADAD). The H» and u methods served as 
the basis for performing the synthesis of the proposed controller. The research [13] considers the p-synthesis 
procedure for developing a robust autopilot. The software-in-the-loop (SIL) verifications applying blade 
element theory (BET) confirms that the autopilot is capable to navigate and land the plane in conditions of 
strong fluctuations in parameters and powerful winds. In studies [14]-[16] Lyapunov functions are used to 
construct robustly stable control systems. The Lyapunov function is constructed in the form of a vector 
function, the anti-gradient of which is set by the components of the velocity vector of the system. Some 
researches [17]-[19] are devoted to the synthesis of robust controllers of aircraft motion parameters using the 
Ha technique. Works [20]-[22] consider the problems of constructing a robust control of the aircraft under 
the action of uncontrolled disturbances, in which the so-called weight functions are introduced. 

Hœ — control theory is widely used in motion control tasks. The modern period of development of 
control theory is characterized by the setting and solution of problems, taking into account the inaccuracy of 
mathematical model of the control object and external disturbances affecting on it. The robust control allows 
to eliminate indicators such as external disturbances and internal parametric uncertainties. The idea of Ho — 
synthesis is to ensure the stability of a closed-loop system not only for a nominal (without model errors) 
object, but also for a "disturbed" object (taking into account the model uncertainties and disturbances 
affecting on the control object) [23]—[25]. 

This paper investigatesH, and Hæ synthesis techniques for the aircraft flight control problem in the 
vertical plane during landing in conditions of undefined disturbances. The application of H,, — control is 
effective for increasing robustness to uncertainties caused by wind perturbations. This document is organized 
as: section 2 describes the principles of stabilization when using H-controls, the synthesis algorithms for 
H, — optimal and H — suboptimal controls, and builds a mathematical model of the aircraft movement in 
the vertical plane with regard to wind disturbances, section 3 presents a mathematical model of a vortex ring- 
shaped wind microburst. The results of Hz and Hœ synthesis methods for the problem of aircraft flight 
control during landing under uncertain perturbations are presented. The efficiency of Hœ — control is 
confirmed by the results of simulation of the synthesized system, and section 4 presents the main conclusions 
of this article. 


2. RESEARCH METHOD 
2.1. Mathematical model of the longitudinal motion of an airplane 

Kinematic and dynamic variables for the equations of motion of the center of mass of the airplane 
are shown in Figure 1. In Figure 1: X is the drag force; x, y is axes of the coordinate system; Y is the lifting 
force; O is the center of mass of the aircraft; V is the airspeed of the aircraft; V, is the aircraft ground speed; 
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wx is the horizontal component of wind speed; wy is the vertical component of wind speed; æ is the angle of 
attack; 0 is the angle of inclination of the trajectory in the air coordinate system. The dynamic equations of 
the aircraft motion in the vertical plane taking into account wind disturbances in projections on the axes of 
the air coordinate system are set by the system of nonlinear differential (1) [26]—-[28]: 


mV = Tcosa — X — mgsin@ — m(w,cosé + Wysind); 
mVÅ = Psina +Y —mgcosé + m(W,sind + W,cos®); 


Jz®z = Mz; 
V = wz; 
h = Vsinð + W, (x, h) 
AT = —(—AT + KypAd;) (1) 
TyB 
A 
y Aircraft axis 


Figure 1. Airplane coordinate system and variables 


The control variables are the thrust T and the angle of attack a, which depend, respectively, on the 
deflection of the engine throttle and the elevator. Here 6, is the engine throttle deflection from the specified 
value. As a result of linearization [29], [30] the nonlinear model of the aircraft (1) is reduced to a linear 
system of ordinary differential equations in increments, which in matrix form has the form (2): 


x = Ax + Bınw + Bonu, (2) 


where x = (AV, A0, Aw,, Ad, Ah, AT)'= the state vector, 
ea, See NR : ; 
w= (w. K Wwe Wwy) = the wind disturbance vector, 
u = (Aôe, A+)" = the control vector. 
The equation for the measured output y in the state space model in the presence of measurement 
noise Ny is written as: 


y= Cyx + Iyny, 
where C, is the measured output matrix and J, is the identity matrix of the corresponding dimension. Thus, 


the mathematical model of the longitudinal motion of the aircraft with considering external wind disturbances 
in the state space model is described by the system (3): 


f* = Ax + Bınw + Boyt, (3) 


y =Cyx + hyny. 


The vector of controlled outputs z;for a linear model of the longitudinal motion of an aircraft taking 
into account wind disturbances in the state space model (3) has the form (4): 


Zı = CX. (4) 
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Consider a vector of controlled outputs Z which is defined as (5): 


=F ]=[) o 


Combining (3), (4) and (5), a system of equations describing the controlled system is obtained: 


x = Ax + By,w + Bonu, 
Zi = Cx, 
Z = 1,uU, 
y =Cyx + yny. 


(6) 


The system of equations describing a standard object in the state space model for an extended vector 
of controlled outputs Z = (z7,u™)" and an extended vector of external inputs W = (w", n3)" has the form 


x = Ax + B W + Bu, 
Z= Cix + DiW + D42u, (7) 
y= Cx + DW + DzU, 


where 
C 0 0 0 
Bı = [Bin 0], By = Ban, Cy = be = Cy, D = p o| Paz = |; |. Pas = [0 1y],D., = 0. 
u 


2.2 Robust H - control 

The algorithms for solving the problems of building H, — optimal and H,,— suboptimal controls are 
considered in this part. Let a finite-dimensional linear controlled and observed object be identified by 
experimental data in the state space model in the form (8) [23], [24]: 


x(t) = Ax(t) + B,w(t) + Bzu(t) 
z(t) = C, x(t) + Dazu (t) (8) 
y(t) = C2x(t) + Dzıw (t) 


where x(t) is the vector of system states; u(t) is the control vector; w(t) is the uncertainty vector, 
characterizing the inaccuracy of the model; y(t) is the vector of measured outputs; z(t) is the vector of 
controlled outputs of the system; A, B1, B2, C1, C2, D11 and Dz, is the constant matrices of corresponding 
dimensions. 

The structural diagram shown in Figure 2 represents the synthesized system. Matrices K(s) and 
G(s) are the transfer matrices of the controller and control object respectively. The matrix G(s) has the 
structure (9): 


A | By B 

— Gy, a — = l = = 
G(s) = Ga Gol |G | 0 Dz (9) 

C | Da 0 


Figure 2. Stabilization principle when using H-control 
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Robust control is included in the lower subsystem loop [G21G22]. Uncontrolled signals passing 
through the upper subsystem [G,,G,2] must be effectively suppressed by it. Let the matrix of transfer 
functions from input w(t) to output z(t) as shown in Figure 2 of the closed-loop system has the form (10): 


= Du D12 -1 C1 
Twz = A + [B,B2]( — K Fie oD K fe] (10) 


The limitation (rationing) Twz is very important. The functional spaces L, (the space bounded by the square 
of the function) and L,, (the space of essentially bounded functions) are considered to describe and define the 
norms T,,, [23]. 

Robust H, and Hœ» — controls are sought in the form of feedbacks u(t) = Ky(t) such that the signal 
norms in the spaces L, and Læ, respectively, equal to ||Twzll2 and ||T,,7||.. are minimal. 


[Tuell = (SI, er€The(-Jo) - TweGo)}do)/? < o (1) 
Twzllo = supë {TwzGo)}<y (12) 
-—O<W<o 
Where ||-|| is the norm in Hardy functional space; W(jw) = w(p)|p=j is the system frequency response; 


tr is the matrix trace; supo{W(jw)} is the maximum singular value of the matrix W (jæ). 
The value ||Ty,||% means the signal energy, and ||Ty,||.. means its intensity in many physical 
applications. Hence, L, is the space of signals of limited energy, and L,, is the space of signals of limited 
intensity. The energy of the error signal under the worst possible perturbation is minimizing while 
minimizing the norm ||T,,||... According to [24], the H, — control equations can be written in the form of an 
optimal observer and an optimal control shaper. 
a = AX(t) + Bzu(t) + L2(C, X(t) E y(t)) (13) 
u(t) = F2 (t) 


Where F, is the matrix of the gain coefficients of the optimal H, — the control feedback; L, is the matrix of 
the optimal feedback gain coefficients by Hz observation. 
Hœ — control has the form (14): 


R(t) = ARC) + ByWryorse(t) + Bout) + ZooLco(Cp%(t) — y(t)) 
u(t) = F(t) (14) 
Wworst(t) = VY 7BLX X(t) 


where Fa is the matrix of the gain coefficients of the optimal Hœ — the control feedback; Z,,L. is the matrix 
of the optimal feedback gain coefficients by Hœ observation. It can be seen from these equations that unlike 
the H, — observer, the Hœ — observer has anobserver-based compensator structure (because of the 
B,Wyyorst(t) component). The main differences in the control structure are the appearance of a new structural 
component B,Wyo7s¢(t) in the Ha — case and the replacement of Lz by ZeoLio- 

H,-optimal control can be constructed in a finite number of operations. However, it is necessary to 
make a reservation that the real program uses iterative procedures to solve the algebraic Riccati equation, 
which makes this statement true if the procedure for solving the Riccati equation considered as a separate 
operation. The Hz-optimal control synthesis algorithm has a linear structure as shown in Figure 3. Unlike the 
H, — case Hœ — suboptimal control (like Hœ — norm) cannot be defined by a finite number of operations and 
requires an iterative procedure. 

The synthesis algorithm of Hœ — control, shown in Figure 4 (in appendix), has a branched structure, 
this is explained by the need to check the condition p(X.Y.) < y? and find y with the required degree of 
accuracy £. If the condition is not met, it is necessary to enter a new value of y larger than the previous one; 
if the accuracy condition |p — po| < £ is not met, where pọ is the spectral radius at the previous value of y 
and p is for current value, it is necessary to enter a new y smaller than the previous one. This algorithm 
implements the input of y at each step manually. Construction of the control is carried out already at the 
selected value of y and the corresponding matrices and X» and Y». It is seen that the synthesis of Hæ — 
control is much more labor-intensive than the synthesis of H, — control also because it is necessary to solve 
two Riccati equations in each cycle of choosing y, while for the H, — case these equations are solved once. 
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Start 


State space matrix input 
A, Bi, B2, Ci, C2 


Decision 
A™X2+X2A-X2B2B2TX2+CiCi7=0 


Decision 
AY2+Y2AT-Y2C2'C2¥2+B1Bi'=0 


Controller design 
F2=-B2™X2; L2=-Y2C2!; 
A=A+BoF2+1L2C>: 
B=-L2; C=F2; D=0 


Figure 3. Flowchart of the H, — optimal control synthesis algorithm 


3. RESULTS AND ANALYSIS 
3.1. Mathematical model of a wind microburst in the form of a vortex ring 

Wind shear is a change in the magnitude and direction of the wind speed during the transition from 
one point in space to another, related to the distance between the points under consideration. The aircraft is 
strongly disturbed during wind shear, therefore, in the take-off and landing modes (due to the low flight 
altitude), a dangerous situation can quickly arise without vigorous and timely intervention of the pilot in the 
aircraft control. Piloting an aircraft in wind shear conditions is challenging due to the need to simultaneously 
control the elevator (balancing the aircraft) and the engine thrust (recovering speed). All of this has led to the 
fact that wind shear at low altitude has become a dangerous phenomenon. Mathematical models of wind 
microbursts of varying complexity have been developed for research and development of automatic or 
manual control systems. 

The generated wind profiles of this model allow simulating atmospheric conditions corresponding to 
some real-life situations, in particular, at Dallas (1985) and New Orleans (1982) airports. According to this 
model [3], [30], the wind microburst area is formed by the flow around a vortex ring located above a flat 
surface. The geometric relations are shown in Figure 5. The mathematical model of a vortex ring-shaped 
wind microburst is described in (15): 


p ==, +r)[F@ — FOI, (15) 


where T is the circulation, R is the radius of the filament of the vortex ring, F(A), F(A) are the complete 
elliptic integrals of the first and second kind, 1,,r, are the largest and smallest distances from the current 
point (x,z,h) to the filament of the vortex ring, Re is the effective radius of the vortex ring core and a 
dimensionless variable. 


n= 


n +T 


The final formulas for wind speeds at the points in space with coordinates (x, h) are as: 
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_ 1,128pP R hp hp R (= hm 
% 2m (((x-x)2-+n24R2)1/2 \r2p Tip ((x-X)2+h2,4R2)1/2 ‘tom Tim f 


_ wef R ( x4 X2 ) R x1 X2 | 
y7 I 24424p2)3/4 \ 3/4 3/4} a 242 4 p2)3/4 \.3/4 3/4 
20 G(x-x)? +h +R?) / rij rp (GO-X)? +h tR?) / Tip Tp 


where 
= _ 358) 2 = y2 2 
Xi =x-X—-R,h, =h—H,r%, = Xf + hb, T2p = x5 + hp: 


To 
X2 =X — X +R, hm =h + H, tym = x2 + h, fom = x3 + hž„ fo = min(Ty,%)),p = 1 — efe. 


The graphs of the vertical w, and horizontal w, components of the wind profile relative to the vortex center 


in the wind microburst zone at a flight altitude of 400 m with parameters (T=45000 m?/s, R=350 m, R,=35 m, 
H=700 m) are shown in Figure 6. 


N 


$ 


Wx, Wy (m/s) 


Time (s) 
Figure 6. Graphs of vertical wy and horizontal w, components of the wind gust profile 
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3.2 Simulation results 

This part considers the specific trajectory of the aircraft glide path [30]-[32]. This trajectory in 
coordinates of altitude and range is a straight line with a given trajectory inclination angle Og; (Ogi= -2.7 
degrees). The task of the control system is to provide a constant air speed Vy = 71.375 m/s and a given height 
h = 400 m when moving along the glide path under the action of wind disturbances, the model of which is 
described above. 

The results of a comparison of the quality of transient processes of closed-loop systems with the 
above described H} — and Hæ — controls constructed using the proposed methodology are presented further. 
During the simulation, the same signal was applied to the input of each closed-loop system, simulating the 
wind disturbance Wacting on the aircraft when it moves in the zone of the wind microburst. Figure 7 shows 
the graphs of the deviation of the airspeed V from the nominal value for two controls. The maximum speed 
deviation when using Hz — control is about 3.49 m/s, and when using Hœ — suboptimal control is about 
1.375 m/s. Hə — suboptimal control provides better quality of transient processes than H, — control 
according to the airspeed deviation. 

Figure 8 shows the graphs of the deviation of the altitude h from the nominal value for two controls. 
Hœ — suboptimal control also provides better quality of transient processes than Hz — control by deviation of 
altitude. Then the maximum height deviation when using H, — control is about 18.75 m, and when using 
Hœ — suboptimal control is about 7.7 m, that is, almost 2.5 times less. This characteristic is very important 
because a sharp loss of altitude in the microburst zone is the main cause of accidents during aircraft landing. 

Thus, it could be concluded that the optimal systems synthesized by the quadratic quality criterion 
are sensitive to the model parameters of the real object and the characteristics of the input influences, i.e are 
not robust [32], [33]. The optimal control has been synthesized by taking into account the need to provide a 
compromise between the minimum possible deviation of the controlled outputs (airspeed and altitude) from 
the nominal values and the power limitations of the control units (engines and elevators). For this purpose, a 
value characterizing the control was introduced into the optimality criterion. H — control solves the 
problem of minimum sensitivity of the closed-loop system for the worst-case external perturbation. The 
energy of the disturbance passing to the output is determined by the Hə — norm of the matrix transfer 
function of the closed-loop system from the external disturbance to the controlled output. The idea of robust 
control synthesis is to provide with one control the stability of a closed-loop system not only for a nominal 
(without model errors) object, but also for a "perturbed" object (taking into account model uncertainties and 
perturbations acting on the control object). 


Time (s) Time (s) 


Figure 7. Velocity deviation V when using Hz — and Figure 8. Altitude deviation h when using H, — and 
H» — control H» — control methods 


4. CONCLUSION 

This paper examines H, and Hœ synthesis methods to reduce the influence of low-altitude wind 
shear on the longitudinal motion of the aircraft in the glide path landing mode. A mathematical model of 
aircraft movement in the vertical plane with considering wind disturbances, and a mathematical model of a 
wind microburst in the form of a vortex ring were obtained. Both controls allow for a significant reduction in 
altitude deviation. However, the Hœ — suboptimal control method provides better quality of transient 
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processes both in altitude and speed than the H, — optimal control. Consequently, the Hœ — control provides 
significantly better suppression of external wind disturbance. Further research by the authors will focus on 
the use of mixed H,/H,, — control, which allows to obtain the average quality of regulation. Despite a much 
more complicated algorithm for its calculation, difficulties in selecting weighting coefficients and the level 
of y, it will probably provide a wide range of transients, each of which in certain casesmay be more useful 
than the results of optimization by a single criterion. Eventually the implementation of mixed H/H» — 
control may create a more comfortable on-board environment than the Hœ — suboptimal system. 


APPENDIX 


State space matrix input 
A, Bi, B2, Ci, C2 and accuracy € 


Decision 
ATXa+XeA-X~(B2B21-y?BiBi1)Xx+C1Ci7=0 
Decision FARE 
AYatYoAT-Yx(C2™C2-y72C17C1) ¥2+BiB17=0 


Calculating the spectral radius 
p = p(Xx, Yo) 
Yes 
Yes 
Controller design 
F..=-B27Xx; La=-YxC2?; Zo=(I-y?7YoXx)"}; ; 


A=At+y?B1Bi'XatB2FotZe2L2C2} 
B=-ZzL2; C=F2; D=0 


Figure 4. Flowchart of the H,, — suboptimal control synthesis algorithm 
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